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ABSTRACT 

We  provide  sufficient  conditions  under  which  two  random  vectors  could 
be  stochastically  compared  using  the  standard  construction.  These  conditons 
are  weaker  than  those  discussed  by  Arjas  and  Lehtonen  (1978)  and  Veinott 
(1965).  Using  these  conditions  we  present  extensions  of  (i)  a  result  of 
Block,  Bueno,  Savits  and  Shaked  (1984)  concerning  the  stochastic  mono¬ 
tonicity  of  independent  and  identically  distributed  random  variables  condi¬ 
tioned  on  their  partial  order  statistics,  and  (ii)  a  theorem  of  Efron  (1965) 

A 

regarding  an  increasing  property  of  Polya  frequency  functions.  Applica¬ 
tions  of  these  extensions  are  also  pointed  out. 

Key  words:  Random  vectors,  stochastic  ordering,  negative  dependence, 

Polya  frequency  function,  component  cannibalization. 


1.  INTRODUCTION 


Let  X  and  V  be  two  real  valued  random  variables  with  survival  functions 
F  and  G  respectively.  Then  X  is  said  to  be  stochastically  larger  than  Y  if: 

(1-1)  F(t)  £  G(t)  ,  t  e  R 

and  is  written  X  Y  [When  equality  holds  in  (1.1)  for  all  values  of  t  we 
write  X  Y.  That  is,  they  are  equal  in  law].  Once  F  and  G  are  given  it 
is  usually  easy  to  verify  (1.1).  The  natural  extension  of  (1.1)  for  finite 
or  infinite  dimensional  random  vectors  is  as  follows  [e.g.  Kamae,  Krengel 
and  O'Brien  (1977)]:  Let  X  =  (X^X-,...)  and  Y  =  (Y^Y-,...)  be  two 
n-component  random  vectors  [n  M  or  n  =  ®]  with  survival  functions  F  and 
G  respectively.  Then  X  is  said  to  be  stochastically  larger  than  Y  if 

(1.2)  P{X  e  A]  £  P { Y  e  A]  for  every  increasing  set  A  e  Rn  , 

and  is  written  X  §st  Y  [When  equality  holds  in  (1.2)  for  all  increasing  sets 
A  e  Rn  we  write  X  S*  Y ] .  In  this  paper  'increasing'  stands  for  'nonde¬ 
creasing'  and  'decreasing'  for  'nonincreasing'.  Unfortunately,  even  with  F 
and  G  explicitly  specified  it  is  usually  not  very  easy  to  verify  (1.2). 

However,  in  light  of 

Lemma  1.1:  X  Y  if  and  only  if  there  exist  two  random  vectors  X 

and  Y  defined  on  a  common  probability  space  such  that  P[X  ^  Y]  =1 

.  A  ,  A 

and  X  -  X  and  Y  _  Y, 


[e.g.  Kamae,  Krengel  and  O'Brien  (1977)],  effort  has  been  made  to  stochas¬ 
tically  compare  random  vectors  by  constructing  them  on  a  common  probability 
space.  In  this  respect  three  alternative  constructions  have  been  used. 

They  are  (i)  standard  construction  [e.g.  Arjas  and  Lehtonen  (1978)],  (ii) 
non-homogeneous  Poisson  process  (NHPP)  construction  [e.g.  Shaked  and 
Shanthikumar  (1984)1,  and  (iii)  total  hazard  construction  [e.g.  Norros 
(1984),  Shaked  and  Shanthikumar  (1985)].  Since  we  will  be  using  the 
standard  construction  we  will  describe  it  here.  Let 

F -j (t)  =  P{X1  >t]  ,  t  e  R 


and 


•  PfX.  >  tlX,  =  xv  X2  =  x2 . X  =  X  }, 


x1,x2,...,Xj_1,  t  £  R;  j  i  2  . 


Note  that  the  dimension  of  the  vector  Xj_^  will  vary  depending  on  where  it 
is  used.  In  Fj(-  f x._1 ) ,  x._1  will  represent  (x1  ,x2, . . .  ,x._1).  We  will  follow 
this  convention  throughout  this  paper.  Define  G1  and  G.(«l^.  1 )  similarly. 

Standard  Construction:  Let  I  { H }  be  the  inverse  function  of  a  survival 

function  H  [that  is,  l{H}(u)  =  inf{t  e  R:  H(t)  >  u],  u  e  (0,1)]  and 
U  =  (U-j/Ug/.-.)  be  an  n-vector  with  independent  components  uni- 

A  A  A 

formly  distributed  in  (0,1).  Construct  X  =  {X^,X2 


, . . .  )  such  that 
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X1  =  l{F1}(U1) 

and  given  X1  =  x1 ,  =  x2,.../Xj._1  =  x._1  , 

X.  =  I  { F.(-  lx.  .)}  (U.),  j  ^  2  . 

j  1  j  — j-i  *  V 

Then  one  has  [e.g.  Arjas  and  Lehtonen  (1978)]. 

Lemma  1.2:  Let  X  be  the  values  obtained  through  the  standard  construc¬ 

tion  . 

Then  X  X. 

Using  the  standard  construction,  Lemma  1.1  and  1.2,  Arjas  and  Lehtonen 
(1978)  [also  see  Veinott  (1965)]  have  obtained  sufficient  condition  on  F  and 
G  that  satisfy  (1.2).  Specifically  they  have 

Lemma  1.3:  Suppose 

IFn(t)  £  G1  (t)  ,  t  e  R 

and 

FjCUx...,)  S  G.dl^.^)  ,  teR  , 

x(  ■>  y.,  i=1 ,2 - j-1;  j  t  2  . 

Then 

(1.4)  X  ^t  Y  • 

One  may  easily  verify  that  if  X  and  V  are  constructed  using  a  common 

A  A 

U  for  both  in  the  standard  construction,  X  £  Y.  Conditions  different  from 
(1.3)  that  imply  (1.4)  have  been  obtained  using  the  NHPP  construction  [see 
Shaked  and  Shanthikumar  (1984)]  and  total  hazard  construction  [see  Norros 


(1984),  Shaked  and  Shanthikumar  (1985)].  In  this  paper  we  provide  con¬ 
ditions  weaker  than  (1.3)  that  imply  (1.4)  and  prove  it  using  the  standard 
construction  [see  Section  2].  Using  these  results  we  obtain  extensions  of 
(i)  a  result  of  Block,  Beuno,  Savits  and  Shaked  (1984)  concening  the  sto¬ 
chastic  monotonicity  of  independent  and  identically  distributed  [i.i.d] 
random  variables  conditioned  on  their  partial  order  statistics,  and  (ii)  a 
theorem  of  Efron  (1965)  regarding  an  increasing  property  of  Polya  fre¬ 
quency  functions  in  Sections  3  and  4,  respectively.  Finally  in  Section  5  we 
establish  the  negative  association  of  i.i.d  random  variables  conditioned  on 
the  partial  order  statistics  and  point  out  a  sample  application  for  results  in 
Section  3. 

The  following  preliminaries  will  be  required  in  Sections  3,4,  and  5. 
Definition  1.1:  A  random  variable  Z  or  its  density  function  h  is  said  to 
be  Polya  frequency  function  of  order  two  [or  log-concave  and  written 
PF^]  if  h(s+x)/h(u+x)  is  decreasing  in  x  e  R  for  all  s  >  u  [Karlin 
(1965)]. 

The  following  closure  property  of  PF^  densities  will  be  needed  later. 
Lemma  1.5:  A  convolution  of  two  PF^  densities  is  PF^. 

Definition  1.2:  Two  random  variables  X  and  Y  with  density  functions  f  and 
g  are  said  to  be  ordered  in  the  sense  of  likelihood  ratio  [and  written 
X  =£pY ]  if  f(t)/g(t)  is  increasing  in  t  e  R  [Karlin  (1965)]. 

Lemma  1.6:  Suppose  X  =£pY.  Then  X  £  Y. 

Definition  1.3:  A  random  vector  T  =  (T^,  T^, . . . ,  Tn)  is  said  to  be  nega¬ 
tively  dependent  through  stochastic  ordering  [NDS]  if 

f(T.,...,T.  -,,T.  T  )l  T.  =  t]  stochastically  decreases  in  t  for 

1  i  -1  i+l  n  i 

all  values  of  i  =  1,2, ...,n. 

Block,  Savits  and  Shaked  (1985)  then  show 
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Lemma  1.7:  If  T  is  NDS,  then 


P{T..  >  t..,...,T  >  t  }  i  U  P{T.  >  t.} 

1  V  '  n  n1  1  i  r 


P{T1  g  S  t  }  S  n  P{T.  g  t.}  . 

1  V  '  n  n  j=1  i  i 

2.  CONDITIONS  FOR  STOCHASTIC  ORDERING  OF  RANDOM  VECTORS 


In  this  section  we  provide  conditons  weaker  than  (1.3)  that  imply 
(1.4)  and  prove  it  using  the  standard  construction.  For  a  given  pair  of 
survival  functions  F  and  G  define 

a-,(t)  =  l{F1}(G1(t))/  t  s  R  and 

l«j.1(ilj-1))KGj.(tly._1))/ 
y1  ,y2'  ‘  * '  ,yj-1 ,t;  e  R  '  J  -  2  ' 

where  a . _ ^  ( y . _ ^  )  and  will  have  different  number  of  components  depending 

on  where  it  is  used.  ct._1(y._1 )  in  F.( . Ia._1  (y._1 ) )  is  (o^  (^  ) ,a2(y2) , .  . .  ,a.  ^ 
(Yj.-]))'  Yj  in  “jCYj)  is  (y1»y2»---/yi)  and  in  cr.(y._1#t)  is  (y1  ,y2,  • . . ,  y-_1 )  ■ 
Then 

Theorem  2.1:  Suppose 


Fq(t)  S  G^(t)  ,  t  e  R 


(2.1)  <  and 


FjCtlofj_1(Yj_1  >)  ^  ^(tl^.^),  yv...,yj_vt  e  R; 


Then  X  §  ,  Y  . 

-  st  - 

Proof:  Constructing  two  random  vectors  X  and  Y  according  to  the  standard 
construction  with  a  common  U  one  sees  that 

A  A  A  A 

(2.2)  X.  =  cr.(Yr  Y2 . Y.)(  j  i  1. 

With  Condition  (2.1)  one  can  easily  verify  that 


v'-'.v's.V.*. -N. 1  A-.  V 
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(2.3) 


a..(t)  £  t,  t  s  R 


=  1  '  VV  y2'**-'yj -v  t£R;  J  = 


From  (2.2)  and  (2.3)  one  sees  that  X  ^  Y  and  the  result  now  follows 
from  Lemmas  1.1  and  1.2. 


Remark  2.1:  From  (2.3)  it  is  immediately  clear  that  (1.3)  implies  (2.1). 

The  popularity  of  using  Lemma  1.3  for  stochastic  comparison  of  random 
vectors  is  its  relative  easiness  to  verify  (1.3).  As  we  will  see  in 
Section  3  there  are  interesting  examples  which  satisfy  (2.1)  but  do  not 
satisfy  (1.3).  In  such  an  example  we  will  also  see  that  it  is  not  hard 
to  verify  (2.1). 

In  some  applications  it  is  possible  to  identify  stronger  conditions  on  F  and 
G  that  imply  (2.1).  Such  a  condition  is  given  in  [see  Remarks  2.2  and 
2.3]: 

Theorem  2.2:  Suppose  there  exist  a  set  of  non-negative  functions 
[b^ ,b^, . . . }  such  that 


(2.4)  <  and 


F1  (t)  ^  G1  (t)  ^  F-j (t+b^ ) ,  t  sR 


F.(tlx.  J  1  G.(tly.  ,)  ^  F.( t  +  b.(x.  .,y.  ,)lx.  -), 

I  -i-1  l  *-1-1  i  i-i-l  *-i-1  -i-1 


t  £  R,  y,  ♦  b,  S  X,  S  y,,  y,  *  S 


xi  s  V|'  1  *  t-2 . M;  i  4  2> 


where  b^  is  a  constant  and  b.  (x^,  y^)  -  b.(x^  . . . /X^ 

yry2, . .  .  ,yM),  i  =  2,3,...  Then  X  £gt  Y. 


Proof;  Observe  that  F^(t)  ^  G^Ct)  i  F^Ct+b^),  t  e  R  implies  t  +  b^  ^ 

a-jCt)  H,  te  R.  Therefore  F2( 1 1 )  §  G^tly.^  =  G2(t  +  b2(x1 , )lx1 ) , 
t  e  R,  +  b^  ^  ^  implies 

(2.5)  F2(tla1(y1)>  S  G^tly.p  ^  F2(t  +  b2(a1(y1  ),y1  )la1  (yn  )),y1  ,teR . 
and  hence 

(2.6)  t  +  b2(a1(y1),y1)  J  o^ty^t)  £  t,  y1 ,  t  e  R  . 

Now  as  an  induction  hypothesis  assume  that 

(2.7)  F.(tlaj_1  (^j_1 ))  ^  G.(tl^j_1 )  ^  Fj(t+b.(a._1(y._1)/^._1)l0(j.1(y..1)) 
and  hence 

(2.8)  t  +  bj(aj_1(y_._1)/yj_1)  £  a{(yM,t)  -  t'  Yy  Y  ■  >Y  t  e  R, 

for  all  i  =  2,3,...,j-1.  Note  that  (2.7)  and  (2.8)  are  valid  for  i  =  2 
[see  (2.5)  and  (2.6)].  Now  from  (2.8)  with  i  =  j-1  and  (2.4)  one 
sees  that  (2.7)  is  satisfied  with  i  =  j.  However,  (2.7)  with  i  =  j 
implies  (2.8)  with  i  =  j.  Then  by  the  induction  hypothesis  one  sees 
that  (2.7)  and  (2.8)  are  true  for  all  values  of  i  =  2,3,...,n.  That  is 
Condition  (2.1)  of  Theorem  2.1  is  satisfied.  The  result  now  follows  by 
Theorem  2.1. 

Remark  2.2 :  As  we  have  noted  in  the  proof  of  Theorem  2.2,  Condition 
(2.4)  implies  (2.1). 

Remark  2.3:  If  we  set  b.  =  <»,  j  Z  1 ,  Condition  (2.4)  reduces  to  Condition 

(1.3).  We  will  use  the  above  result  to  obtain  an  extension  of  a  Theorem  of 


Efron  (1965). 
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3.  STOCHASTIC  MONOTON ICITY  OF  i.i.d  RANDOM  VARIABLES 


CONDITIONED  ON  ORDER  STATISTICS 


Let  Z.  =  (Z.j /Zg, . . .  ,Z.)  be  a  random  vector  of  j  i.i.d  random  variables 
with  survival  function  H  and  density  h.  Now  let 
^.(t)  =  #{i  :Zj  >  t,  i  =  1,2,. 

{ N j ( t) }  is  assumed  to  be  right  continuous  with  left  hand  limits.  For  a  given 

r  £  1 ,  z  =  (z1  ,z2, . . .  ,zr),  i  =  ,£2, . . .  ,lr)  and  m  =  (m^  m2,...,mr) 

such  that  z^  <  z £  <• .  -<zr  and  0  S  SLV  £  m^  l..i  1^  =  ™2  *  -  mi  =  J'  de" 

fine  Z.  „  =  (Z.IN.(z.)  =  £.,  N.(z.-)=m.,  i  =  1,2, ...,r} 

-j,z,£,m  l-j  j  1  1 '  J  *  1 


Let  T  =  Z.  „  .  Then  simple  calculation  shows  that 

-j,z,.2,m  K 


(3.1)  P{T1  >  t}  =  .  {m.  +  UM  -  mj)H2  .(t)} 


where 


Hz,i^  = 


H(t)  -  H(z.) 
H(zh1)  -H(Zj) 


z.  1  <  t  <  z., 


i  =  1,2, ... ,r  +  1 

and  Zq  =  zr+1  =  +  °°.  Since  .(t)  increases  as  Zj_^  and  z ■  increases 
for  Zj_i  <  t  <  Zj,  it  is  not  hard  to  see 

Lemma  3.1:  P >  t}  given  by  (3.1)  is  increasing  in  (z,  £,  m). 

Now  consider  P{T2  >  tlT^  =  t^}.  A  routine  calculation  shows  that 


(3.2)  P{T„  >  tlT,  =  tj  = 


(m/  +  (Si;  ,  -  m/)  H_  ,(t)}. 
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where 


ii 

'<W  j 

(3.3) 

|  and 

i 

t 

1 

i  m*j  *  ", 

■'(z.st^  '  i  =  1'2 . r-  : 

Define 

I'  =  (Tv  T^,.. 

.  T.'_1 )  =  Z.  ,  .,  ,.  Then  comparing  the  right  ; 

j  ” '  j  *  *  /  f  i 

hand  side  of  (3.2)  and  that  of  (3.1)  one  sees  that  P{T2  >  tlT^  =  t^}  = 
P{T'  >  t},  t  e  R.  Now  considerng  P{T3  >  tlT.j  =  t2}  and  continuing  this 
one  can  establish 


Lemma  3.2: 


st 


«V  X3 - T.X-I-,  =t,l  = 


, ,  where  £'and  m' 


are  given  in  (3.3), 


Remark  3.1:  Note  that  the  event  {Nj(zp  =  N^(Zj-)  =  m.,  i  =  1,2,...,r} 

relays  information  on  the  order  statistics  z^j  ^  Z[2]  =  -  -  =z [ j ]  of  Z.y 
For  example  consider  the  following  event  [1  £  k1  <  k2  <...<kr  S  j]: 


(3.4)  A  =  { N . ( 2j )  =  j  -  kj7  N^Z;-)  =  j  -k.  +  1,  i  =  1,2, ...,r}, 

for  some  1  $  r  %  j . 

Then  clearly 

A  E  ^ Z [ k^ ]  =  ZV  Z[kz]  =  Z2'""Z[kr]  =  zr}- 

Block,  Bueno,  Savits  and  Shaked  (1984)  consider  the  random  vector 
conditioned  on  event  A  as  given  in  (3.4)  and  study  its  stochastic 
monotonicity  and  NDS  properties.  Alternatively  consider 

(3.5)  A  =  {N.(Zj)  =  j-k.+l  =  Nj(z.-),  i  =  1,2,. ..,r}. 

Then 

A  =  <  Zj  <  i  =1,2,  ...,r},  where  we  set  Z^  = 

Now  we  present  the  main  theorem  of  this  section. 
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Theorem  3.3:  Let  Z  .  _  be  as  defined  in  this  section  with  j  =  n.  Then 
-  — n,z,£,m 

Z  ,  o  m  is  stochastically  increasing  in  (z,  £,  m)[that  is  (w,u,v)  t 

(c,a,b)  =>  Z  £  .  Z  .]. 

--  -  -n,w,u,v  st  -n,c,a,  b J 

Proof:  Let  F  and  G  be  the  survival  functions  of  X  =  Z  and 

-  -  -n,w,u,v 

Y  =  Z  .  ,  respectively  for  some  (w,u,v)^  (c,a,b).  From  Lemma 

n  /  C  /  a  r  D  * 

3.1  one  has 


(3.6) 


F^t)  £  G^(t),  t  e  R 


Now  consider  a^(t)  as  defined  in  Section  2.  From  (3.6)  it  is  clear 
that  a^(t)  g  t.  A  careful  study  of  F^  and  G^  [given  in  (3.1)  with 
appropriate  values  for  z,  £  and  m]  will  show  that  for  fixed  i,  Cj  and 


(3.7) 


If  u.  =  a.  then  {t:t>c-}  =  (t:a-(t)  >  w.}, 

i  i  i  1  1  i 

If  Uj  >  af  then  [t:t  >  c.}  C  {t:o^(t)  >  w.}  , 


if  v.  =  b.  then  { t : tSc- }  =  [t:a-(t)  £  w.},  and 
i  i  ‘  r  1  r 

if  Vj  >  bj  then  {t:t  £  Cj}  C  [tia^t)  S  w.}. 


Therefore  if  we  set 

Ui  '  Ui  ' ( w . <  a^(yn))  '  vi  “  vi  ’(w.S  a^(y1  )) 


ai  ~  ai  *(c  <  y  )  '  bi  bi  '(c.^  Y*  )  ' 
from  (3.7)  one  sees  that  for  all  y^  e  R,  (u',  v')  £  (a',  b').  There¬ 
fore  from  Lemmas  3.1  and  3.2  one  has 


R 


*  2'  3' * ’  ’ '  n  1 


T  v  y  1  J 


-n-1 , w,u ' ,v' 


(VV3 . YnIY1  =  vn}  51  ln^,c,a’,b' 


and 

FgCtla^Cy^))  i  G2(tly1 )  ,  y1 ,  t  e  R. 

Therefore  a  continued  application  of  the  above  analysis  will  result  in 
-1(^-1))  =  ^(t!^),  y1,y2,.../yj.1/  t  e  R  ;  j  £  2  . 

That  is  (F.)  and  (G.)  satisfy  Condition  (2.1)  of  Theorem  2.1  and  there¬ 
fore  X  §st  Y  . 


Now  combining  (3.4)  with  the  above  result  one  obtains 
Corollary  3.4:  For  some  H  r  S  n,  z1  <  <...<zr,  and  H  <  kg  <... 

<  kr  £  n  let  T  =  *  V  z[k,]  =  *2 . z[krl  *  *r>'  Then  - 

is  (i)  stochastically  increasing  in  z  and  (ii)  stochastically  decreasing  in 
k. 


Block,  Beuno,  Savits  and  Shaked  (1984)  proved  (i)  [see  the  Corollary 
to  Lemma  4.2  there].  They  also  used  this  result  to  show  some  negative 
dependence  properties  of  T.  In  Section  5  we  will  show  the  negative 
dependence  property  of  Z  0  .  Now  combining  (3.5)  with  Theorem 

3.3  one  obtains  [as  a  special  case]. 
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Corollary  3.5:  For  some  1  S  k  £  n  and  z  e  R  let  T  =  { Zn  I Z  ^  ^ 


<  z  < 


Z^j}.  Then  T  is  stochastically  increasing  in  z  and  is  decreasing  in  k. 


We  will  use  the  above  result  to  provide  a  decision  rule  for  selecting 
components  during  cannibalization  [see  Section  5]. 


4.  INCREASING  PROPERTY  OF  POLYA  FREQUENCY  FUNCTIONS 

Let  Z  =  (Z^,  Z2/...,Zn)  be  n  mutually  independent  PF2  random  vari¬ 
ables  with  density  functions  (h^,  h2#...,h  ).  Let  Zs  =  [ZIZ  2=sl  /  where 
V  =  (1,1  ,...,1)  is  of  appropriate  dimension. 

Theorem  4.1:  For  some  s  i  u,  let  F  and  G  be  the  survival  functions  of  Zs 
and  Zu  respectively.  Then  (i)  F  and  G  satisfy  Condition  (2.4)  of  Theorem  2.2 
and  (ii)  Zs  *st  Zu. 

Proof:  From  Theorem  2.2  and  (i),  (ii)  follows.  So  consider  (i).  Let  be 
the  density  function  of  the  sum  Zj  +  Zj+1  +...+  Zn<  Then  if  (f.)  and 
(g.)  are  the  densities  of  (F.)  and  (G^)  respectively, 


f,(t)  1 

1  q2(s_t)  J 

1  ^(u) 

g-,(t) 

1  q2(u-t)  i 

[  q1(s) 

Note  that  s-t  £  u-t  and  since  q.  is  PF2  [Lemma  1.5],  q2(s-t)/q2(u-t)  and 
hence  f ^ ( t )/g^ ( t)  increases  in  t.  That  is  and  from  Lemma  1.6, 

has 

F^(t)  g  G^(t)  ,  t  e  R  .  Now  consider 


f^t  +(s-u)) 

g^t) 


h^(t  +  (s-u)) 

h^t) 


q^u) 

q^s) 


one 


Since  is  PFg,  h1(t-(s-u))/h1(t)  and  hence  f1(t-(s-u))/g1(t)  decreases  in 
t.  Therefore  as  before  V1  2£r  X^Cs-u)  and 
P 1  (t)  £  G-j(t)  2  F^t  +  bp  , 
where 

b^  =  s-u  . 

Now  let 


j-1 


bj(*j-V*j-n)  =  s-u-  Z^  (x.  -  y.)  ,  j  =  2,3, . . . ,n. 


From  (4.1)  one  has  y1  +  b^  2  a^Cyp  2  y^ .  Hence  with  restriction  y1  +  b1 


2  y^  one  sees  that 


b2(x1,y1)  2  0  . 


Now  as  an  induction  hypothesis  assume  that 
( A1 )  b.(x._1  )  2  0  ,  j  S  £  and 

(A2)  Condition  (2.4)  holds  for  all  values  of  j  up  to  £-1. 

Note  that  (A1)  and  (A2)  is  true  for  £  =  2  [see  (4.1)  and  (4.2)].  Consider 


£-1 


£-1 


I  qo+ its-t-  *  x>.) 


£+1 


i=1 


£-1 


g£(tl^£-1))  \  q£+1^u"t_  yi) 


q£(u"  1  yi) 
i=1 


£-1 


q£(s‘  1  V 

*  i=1  1 


Since  by  Assumption  (A1),  b£(x£-1  )  2  0,  one  has 


£-1  £-1 
s-t-  Z  x  >  u  -  t  -  Z  y. 
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Therefore  as  before 


*X£,X1  =  x1,...,X£_1  -  x£_1}  ^£p  {Y£IY1  yv...,Y£_i  -  y£_-,}  , 
Yj  +  t>j(xj_1,  yM)  S  Xj  S  yj,  i  =  1,2,. ..,£-1,  and 

^£***-£-1*  “  ^£(tl*£-1^  • 

Now  consider 

f£(t  *  b£(S£-1'W*£-1>  [  W51.1''!h,))  (  q£(u‘J  ' 


£-1 

9t(t,^-1>  '  h£(t)  '  '  q£<S_  it1Xi) 

Since  by  Assumption  (A1)  b£(x£_1/y£_1 )  S  0,  the  log  concavity  of  h£ 

implies 

(4.3)  F£(t,x£_1)  =.  G£(tly£_1)  5  F£(t  +  b£(x£1  ,y£-1  )lx£_1 )  , 

for  all  y1  +  bq  £  x1  >y1 ,  y.  +  b.(xM , ^  x.  S  y.  ,  i  =  1,2,.. .,£-1. 
With  (4.3)  one  has 

(4.4)  y,,  *  8  «e(^>  8  V£  ■ 

Therefore  with  the  restriction  y^  +  b^  S  x^  ^  y^ ,  y;  +  b.(x._.j,  y^)  - 
Xj  i  y.,  i  =  1,2,..., £-1  one  obtains 

(4.5)  b2+1(x£,y£)  =  b£(x£_1  ,Y£.-|)  -  (x£  -  y£)  £  0. 

From  (4.3)  and  (4.5)  and  the  induction  hypothesis  one  sees  that 
F  and  G  satisfy  Condition  (2.4). 

Remark  4.1:  Efron  (1965)  using  an  alternative  proof  established  (ii) 
of  Theorem  4.1.  The  discrete  analogue  of  this  result  has  proved 
to  be  very  useful  in  the  analysis  of  queueing  networks  [see 
Shanthikumar  and  Yao  (1985a),  (1985b)].  We  will  next  provide  an 
extension  [the  main  result  of  this  section]  of  Theorem  4.1.  This 
extension  has  useful  applications  in  the  analysis  of  queueing 
network  and  these  applications  will  be  discussed  elsewhere. 


Theorem  4.2:  Let  X  =  (X.,,  X£,...,Xn)  and  Y  =  (Y.,  ,  Y^...,  Y  ) 

be  two  independent  random  vectors  of  mutually  independent  elements 

A  A  A  A  A  A 

with  PFg  densities  (  fy  f2 , . . . ,  f  n )  and  (g^,  g2,...,gn)  respectively. 
For  some  s  and  u  define 

A  A  A  A 

X  =  {XIX  1  =  s}  and  Y  =  { VI V  1  =  u}  . 

Suppose  for  some  0  i  aj  <  » 

%Yi  Vi  -  *i  •  1  *  1'2 . "■ 

n 

Then  for  s£u  +  I  a.  ,  X  £  Y  . 

i=1  1  sl 

Proof:  Suppose  we  can  show  that 

<4'6>  fXV  Y2 .  VX,  ...-Yn.v)  ist 

A  A  A  A  A  A 

{YvY2,...,YnIY1  +  Y2  +  ...  +  Yn  =  u}  ,  v  £  u  +  a.,  . 

A  A  , 

Then  applying  (4.6)  with  X2  replacing  Y2  one  sees  that 


A  A 


A  A 


{Xr  X2,  Y3,...,Yn)X1  ^ 


+  X2  +  Y3  +  . . .  +  Yn  =  w} 


A  AS 


A  /V 


V  V3 . W 


AS  <A 


+  Y-  +  Y-  +  . .  .  +  Y  =  v},  w  £  v  +  a 
2  i  r> 


n 


AS  AS 


=st*Y1'  Y2'  Y3,,--,YnIY1^ 


AS  AS 


+  Y2*Y3  ♦...+  Yn.u), 


w  £  v  +  a. 


£  u  +  a1  +  a2. 


Continuing  this  way  one  can  obtain  the  desired  result.  Hence  all  we  need 
is  to  establish  (4.6)  for  v  =  u  +  a^.  So  without  loss  of  generality  let  the 
left  hand  side  of  (4.6)  be  X  and  the  right  hand  side  be  Y,  with  survival 
functions  F  and  G,  and  densities  f  and  g,  respectively.  Then. 
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gi(t)  )j| 


q2(v_t)  )  (  9i  * 

q2(u_t))  f  fl  *  q2(V) 


where  q_  is  the  PF?  density  of  I  Y.  and  *  stands  for  convolution. 

d  d  i=2  1 

As  before  one  sees  that  f1(t)/g1(t)  is  increasing  in  t  arid  hence 
F-j  (t)  £  G1  (t)  ,  t  e  R  . 


Consider 


f.,(t  +  a1) 
g-,(t) 


f.,(t  +  a1 ) 

g-,(t) 


g1  *  q2(u) 

— X - 

f!  *  <?2(v) 


The  above  expression  is  decreasing  in  t  and  hence 
F1  (t)  >  G^t)  £  F1  (t  +  ai). 

So  F^  and  G1  satisfy  (2.4)  with  b1  =.  a1  £  0.  Now  note  that 


{xr  x2,.../xnix1  =  = 


{Y2,  Y3,  .  .  •  /  YnIY2  +  Y3  +  .  . 

v  -  x.,}  =  YV'X1 


+  Y  = 
n 


and  similarly 
^Y2'  Y3' 


<Y2'  y3'"'  VnIY1  *  V,}  = 

<Y2'  Y3 . YnIY2  *  Y3  *  ■■■  *  Yn  = 

u-y,I  s  Yu-V,  . 

Then  for  y^  +  b^  =  x-|  =  ,  v  -  x^  2  u-y^,  and  therefore  from  Theorem  4.1 

one  sees  that  F.  and  G^,  j  £  2  also  satisfy  Condition  (2.4)  of  Theorem  2.2. 

Hence  XL  Y  . 

-  st  - 

A  .  A 

Remark  4.2:  Obviously  when  Xj  -  Y.,  one  can  set  a  j  =  0 ,  i  =  1,2,...,n. 

Then  Theorem  4.2  specializes  to  Theorem  4.1.  Now  suppose  instead  of 

-A  A  A 

conditioning  X  1  =  s  and  Y  1  =  ii,  we  condition  them  to  be  X  T_  -  Z^ 

A  . 

and  Y  1  -  Z2,  where  and  Z2  are  two  random  variables  with  distri- 


WWW 


••**•*.  V*  .’V*  w'-  »'■  .*•  *  *  *  *•■  v.Cvvv 

»\  <  3  •.  •*.  ••  f  .  «.  •  -  *.*.•-*  .«*.^**  I".  V  •  w*.  s.  A. .  V  aU.  AiV  •-  ‘  V  *-  -  ■ 


* 


butions  and  H 2,  respectively.  That  is, 

A  A 

P{X  e  A}  =  /  P{X  e  AIY  1  =  s}  dH^s)  and 

A 

P{Y  e  A}  h  /  P{ Y  e  A  I  Y  1  =  u}  dH2(u)  ,  A  £  Rn  . 

n  n 

Suppose  Z.  5  t  Z~  +  I  a.  .  Then  defining  Z'  =  Z.  -  I  a.  with 
i  st  ^  j=i  i  it  j=1  i 

distribution  function  H.j  one  obtains 

P{X  e  A}  =  f  <£%  (A,  u)  dH^(u)  and 

P{Y  £  A}  =  f  ^(A,u)dH2(u)  ,  A  £  Rn  , 

n 

where  <yj  (A,u)  =  P{X  £  AIX  2  =  u  +  I  a.} 

i=1  1 

^(A,u)  =  P{Y  £  AIY  1  =  u}. 

From  Theorem  4.2  one  has  3^(A,u)  £'^(A,u)  for  all  increasing  upper  sets 
A,  and  from  Theorem  4.1  one  has  ^;(A,u)  and  /'(A,u)  increasing  in  u. 

d 

n 

Then  it  is  easily  seen  that  since  ^st  Z2  +  I  a.f  [that  is  H^(u)  %  H2(u)], 

one  has  P[X  £  A}  S  P[Y  £  A}  for  all  increasing  sets  A  £  Rn.  Therefore  X 

£  «.  Y. 
st  - 

Using  analysis  similar  to  that  in  this  section  it  can  be  shown  that  all 
the  results  presented  in  this  section  holds  true  for  discrete  random  vari¬ 
ables.  A  special  case  of  the  above  extension  has  been  used  by  Shanthi- 
kumar  and  Yao  (1985b)  in  the  analysis  of  closed  queueing  networks. 


5.  APPLICATIONS 


We  will  first  apply  Theorem  3.3  to  show  that  Z  ,  defined  in 
7  -n,z,£,m 

Section  3  satisfies  the  NDS  property  of  Block,  Savits  and  Shaked  (1985) 
[see  Definition  1.51. 


Theorem  5.1:  T  =  Z  „  satisfy  the  NDS  propetry. 

- -  -0,2,2,131  1  r  r  i 

Proof:  From  Lemma  3.2  it  can  be  seen  that 

(5.1)  {(Tr  T2,.../Tj_1,  Ti  +  1/- /Tn)|Tj  =  S1  -  -n-1  ,z,2s,ms'  where 
2®  =  2.  -  l^z  <  and  ms  =  nr  -  l^2  ^  i  =  1,2,...,n.  Therefore 

s  >  t  =  (2s,  ms)  ^  (2*,  mt).  Then  from  (5.1)  and  Theorem  3.3  one 
has  the  NDS  property  for  T. 

Next  we  point  out  an  application  for  Corollary  3.5  in  component  cannibali¬ 
zation.  Consider  a  collection  of  heat  sources  each  cooled  by  its  own  cooling 
system  consisting  of  a  set  of  n  identical  pumps  and  a  circulation  system 
[composed  of  radiators,  pipes  etc.].  The  operation  of  the  heat  source  is 
continued  unless  either  the  heat  source  or  the  cooling  system  fails.  A 
cooling  system  failure  may  occur  either  because  all  the  n  pumps  have  failed 
or  because  the  circulation  system  has  failed.  The  circulation  system  fails 
mainly  due  to  the  structural  damage  caused  by  fewer  number  of  pumps 
working.  One  may  safely  assume  that  the  damage  accumulation  rate  increases 
as  the  number  of  working  pumps  decreases. 

After  t  time  units  of  operation  all  the  cooling  systems  are  replaced  by 
spare  cooling  systems.  This  t  time  units  may  represent  the  ‘high  reliable 
operatng  time1  for  the  cooling  system.  These  pulled  out  cooling  systems 
may  however  be  used  elsewhere  for  less  critical  use.  Since  some  of  the 
pumps  from  these  pulled  out  cooling  systems  may  have  already  failed  one 
could  pool  all  such  cooling  systems  and  assemble  fewer  number  of  systems 
but  all  with  n  operating  pumps.  In  such  a  case  one  has  to  choose  the  best 
circulation  systems  for  these  re-assembly.  Since  it  is  not  economically 
feasible  to  test  the  circulation  system  for  its  structural  damage  one  may  use 
the  following  result  to  choose  the  desired  circulation  systems. 

Let  C.  be  the  circulation  system  belonging  to  the  i-th  cooling  system. 


Assume  that  we  have  also  observed  that  kj  of  the  n  pumps  of  the  i-th 
cooling  system  were  not  operative  at  time  t.  Then  if  Dj  is  the  damage 
accumulated  in  Cj,  one  has  from  Corollary  3.5 

(51>  °n(1)  Sst  D„(2)  Ss,  •••  4  D*(m>  ' 

where  n  is  a  permutation  of  {1,2,...,m}  such  that 

(5'2)  k7i(1)  =  kn(2)  -  ••  ~  kn(m)  • 

One  may  now  use  (5.1)  and  (5.2)  to  choose  the  stochastically  better 
circulation  systems  for  re-assembly. 
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